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Abstract. The electro-diffusion model, which arises in electrohydrodynam- 
ics, is a coupling between the Nernst-Planck-Poisson system and the incom- 
pressible Navier-Stokes equations. For the generally smooth doping profile, the 
quasineutral limit (zero-Debye-length limit) is justified rigorously in Sobolev 
norm uniformly in time. The proof is based on the elaborate energy analysis 
and the key point is to establish the uniform estimates with respect to the 
scaled Debye length. 



1. Introduction and Main Results 

In this paper we consider a model describing ionic concentrations, electric poten- 
tial, and velocity field in an electrolytic solution. This model is a coupling between 
the Nernst-Planck-Poisson system and the Navier-Stokes equations [18, 17, 2, 14]. 
The (rescaled) system takes the form 

= div(Vn^ - n^V$^ - n^w^), (1.1) 

= div(V/ + p^V$^ - p^v^), (1.2) 

X^A<i>^ ^n^ ^p^ - D{x), (1.3) 

Vt + ■ Vv^ + Vtt^ - fiAv^ = {n^ - /)V«>^, (1.4) 

divi;^=0 (1.5) 

with initial data 

n^{x,0) ^n^{x), p^{x,0) =p^{x), v^{x,0) = v^{x), a; G T^ (1.6) 

where T"^ is the periodic domain in M.^ , and p^ denote the negative and positive 
charges respectively, the electric field, the velocity of the electrolyte, and tt^ 
the fluid pressure. The parameter A > denotes the scaled Debye length and /i > 
the dynamic viscosity. D{x) is a given function and models the doping profile. 

Usually in electrolytes the Debye length is much smaller compared the others 
quantities, and the electrolytes is almost electrically neutral. Under the assump- 
tion of space charge neutrality, i.e. A = 0, we formally arrive at the following 
quasineutral Nernst-Planck-Navier-Stokes system 

nt = dYv{yn + n£ — nv), (1-7) 
Pt = div{Vp - p£ - pv), (1.8) 
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n — p — D{x) = 0, 

vt + v ■ Vu + Vtt — /uAt; = — (n — p)f , 
divu = 0, 



(1.9) 
(1.10) 
(1.11) 



where we assume that the Umits n, j/ p, v, — V^"^ = ^ £ exist 

as A ^ 0+. 

The purpose of this paper is to justify the above hmit rigorously for sufficiently 
smooth solutions to the system (1.1)-(1.5). 

Since the incompressible Navier-Stokcs equations (1.4)-(1.5) are involved in the 
system (1.1)-(1.5), it is well known that whether the global classical solution for 
general initial data exists or not is open for three spatial dimensional case and 
only local classic solution is available. For example, in [14], Jerome studied the 
Cauchy problem of the system (1.1)-(1.5) and established the local existence of 
unique smooth solution for smooth initial data. The local existence of unique 
smooth solution to the incompressible Navier-Stokes equations can be obtained by 
standard method, see [12, 19]. 

The local existence of unique smooth solution to the limiting system (1.7)-(1.11) 
with initial smooth data 

n{x,t = 0) = no{x), p{x,t = 0) = po{x), v{x,t = 0) = vo{x) (1-12) 

can be obtained by the similar arguments to those stated in [14]. Since we are 
interested in the quasineutral limit of the system (1.1)-(1.5), we omit the detail 
here. 

In this paper we assume that the doping profile is a smooth (sign-changing) 
function and the initial data nQ{x),pQ{x) and Vq{x) are smooth functions satisfying 



The main result of this paper can be stated as follows: 

Theorem 1.1. Let [n^ ^p^ ,v^) , = — be the unique local sm,ooth so- 
lution to the system (1.1)-(1.5) with initial data (1.6) on x [0,T*) for some 
< T* < 00. Let {n,p, £, v), £ = — V$ he the unique smooth solution to the limiting 
system (1.7)-(1.11) with initial data (1.12) on T'^ x [O.To) for som,e < Tq < +oo 
satisfying n + p > kq > 0, where kq is a positive constant. Suppose that initial data 
satisfy (1.13) and 

n^(,r) = no{x), p^{x) = pa{x) + \^<Xw£{t = 0), v^{x) = va{x). (1.14) 

Then, for any T G (0, min{To, T^}), there exist positive constants K and Aq, Aq <C 1, 
such that, for any A e (0, Aq), 



for any a € (0, 2), independent of A. Here = — n, p^ = p^ — p, E^ = E^ — £, 
and = — v. 




(1.13) 



sup 

0<t<T 



{\\{h\p\E\v^){t)\\H^+\\{n^,pi,i^m\U^ 
+ X\\E\t)\\H^ + X\\E^{t)\\„r } < KX^-^/^ 



(1.15) 
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Remark 1.1. In this paper wc deal with the three spatial dimensional case, if the 
problem (1.1)-(1.5) is considered in two dimensional space, both the problem (1.1)- 
(1.5) and the limiting problem (1.7)-(1.11) enjoy global smooth solutions, thus we 
can obtain a similar result to that stated in Theorem 1.1 (in fact much easier). 

Remark 1.2. If the assumption (1.14) does not hold, we need to consider the initial 
layers. On the other hand, if we consider the system (1.1)-(1.5) on the smooth 
bounded domain in M.^, the boundary layers may appear. These issues will be 
studied in the future. 

The main difficulty in dealing with the quasineutral limits is the oscillatory 
behavior of the electric field (the Poisson equation becomes an algebraic equation 
in the limit). Usually it is difficult to obtain uniform estimates on the electric 
field with respect to the Debye length A due to a possible vacuum set of density. 
To overcome this difficulty, we introduce the following A-weighted Lyapunov-type 
functionals 

+ X^\\{E^ ,d\vE^ ,VdwE^ ,dwE^)\\^ + ||(£;^, div£;^)|p (1.16) 

and 

G^{t) ^ \\{A~zt,AiuEldwE^)\\l.+\^\\yAWE^\\\ (1.17) 

where z^ = -\- p^, = — n,p^ = — p,v'^ = v'^ — v, E^ ~ E^ — £ and 
{n^,p^,E^,v^) denotes the difference between the solution to the system (1.1)- 
(1.5) and the solution to the limiting system (1.7)-(1.11), see Section 2 below for 
details. By a careful energy method, we can prove the following entropy production 
integration inequality 

r^{t)+ f G^{s)ds<Kr^{t = 0)+KXi + K{T^{t)y +K f T^{s)G^(s)ds 
Jo Jo 

+ K I [T^{s) + (r^(s))']ds, t > (1.18) 

for some positive constants q, r, K and I, independent of A, which implies our desired 
convergence result by the assumption of small initial data r'^(O). 

Remark 1.3. The inequality (1.18) is a generalized Gronwall's type with an extra 
integration term where the integrand function is the production of the entropy and 
the entropy-dissipation. Hence (1.18) is called as the entropy production integration 
inequality. 

Remark 1.4. The A-weighted Lyapunov-type functional (1.16) and (1.17) is moti- 
vated by [10, 21], where the quasineutral limit of drift-diffusion-Poisson model for 
semiconductor was studied. However, in our case the incompressible Navier-Stokes 
equations are involved and the more refined energy analysis is needed. We be- 
lieve that those A-weighted Lyapunov-type energy functionals can also be used to 
deal with the quasineutral limit problem of other mathematical models involving in 
Navier-Stokes equations, for example, the mathematical model for the deformation 
of electrolyte droplets: 

p{ut + u ■ Vu) + n = fAu +{n- p)W - V • (Vcj) V(/)), 
V-u = 0, 
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Tit + u ■ Vn 
Pt + u- Vp 

V ■ (Avy) 

(j)t + u-V(t) 



V • (DnVn - /x„nVy + MvNcj)), 

V • (£)pVp - iJippVV + MpV4>), 
n-p, 



^{A<f>-rj-^W'{<f>)), 



where 7, i^, 77, 1?„, fj.n,fj.p and M are positive constants, see [16] for the detailed 
description on this model. 

We point out that the quasineutral hmit is a well-known challenging and physi- 
cally complex modeling problem for fluid dynamic models and for kinetic models of 
semiconductors and plasmas and other fields. In both cases, there only exist partial 
results. For time-dependent transport models, the limit A — > has be performed for 
the Vlasov-Poisson system by Brenier [1] and Masmoudi [13], and for the Vlasov- 
Poisson-Fokker-Planck system by Hsiao et al [8, 9], respectively. For the fluid dy- 
namic model, the drift-diffusion-Poisson system is investigate by Gasser et al [6, 7] 
and Jiingel and Peng [15], and for the Euler-Poisson system by Cordier and Gre- 
nier [3] and Wang [20]. Recently, Wang et al [21, 10, 23] extends some results cited 
above for the general doping profiles, the main idea is to control the strong nonlinear 
oscillations caused by small Debye length by the interaction of the physically moti- 
vated entropy and the entropy dissipation. For the Navier-Stokes-Poisson system, 
Wang [20, 22] obtained the convergence of the Navier-Stokes-Poisson system to the 
incompressible Euler equations. Ju et al [11] obtained the convergence of weak solu- 
tions of the Navier-Stokes-Poisson system to the strong solutions of incompressible 
Navicr-Stokes equations. Donatelli and Marcati [4] studied the quasineutral-type 
limit for the Navier-Stokes-Poisson system with large initial data in the whole space 
R'^ through the coupling of the zero-Debye-length limit and the low Mach number 
limit. 

We mention that there are a few other mathematical results on the system (1.1)- 
(1.5). Jerome [14] obtained the inviscid limit {fj, 0) of the system (1.1)-(1.5). 
Cimatti and Fragala [2] obtained the unique weak solution to the system (1.1)-(1.5) 
with Neumann boundary condition and the asymptotic behavior of solution when it 
is a small perturbation of the trivial solution for the stationary problem. Feireisl [5] 
studied the system (1.1)-(1.5) in periodic case without the diffusion terms in the 
first two equations and obtained the existence of weak solution. 

Before ending this introduction, we give some notations. We denote || • || the 
standard norm with respect to x, the standard Sobolev space W'^''^, and 
II ■ W^k the corresponding norm. The notation ||(^i,A2,''' i^n)|P means the 
summation of ||^i|p,i = I,-- - ,n, and it also apphes to other norms. We use 
Cj, 6i, e. Kg, Ki, and K to denote the constants which are independent of A and 
may be changed from line to line. We also omit in integral spatial domain for 
convenience. In Section 2, we give some basic energy estimates of the error system, 
and the proof of Theorem 1.1 is given in Section 3. 



In this section we obtain some energy estimates needed to prove our result. To 
this end, we first derive the error system from the original system (1.1)-(1.5) and 
the hmiting system (1.7)-(1.11) as follows. Setting = — n,p^ = p^ — p,v^ = 

-v,'7r^ = 7r^ -Tr,E^ =E^ -£ with = -Vl>^, E^ = -V$^, £ = -V$ and 



2. The energy estimates 



QUASINEUTRAL LIMIT OF THE ELECTRO-DIFFUSION MODEL 5 

|,A _ _ ^ using the system (1.1)-(1.5) and the system (1.7)-(1.11), we obtain 

ht = div(Vn^ + nE^ + fi^{E^ + S) - n^{v^ + v) - nii^), (2.1) 

= div(Vp^ - pE^ - p^{E^ P^{v^ + v)- pv^), (2.2) 

- X^divE^ = fi^ -p^ + A^divf, (2.3) 
v^ + ■ Vv^ + V ■ Vw^ ■Vv + Vtt^ - iiAv^ 

= -{n^-p^){E^+S)-{n-p)E\ (2.4) 

di-vv^ = 0. (2.5) 

Set Z = n + p, then (1.7)-(1.11) is reduced to 

Zt = dw{yZ + D8 - Zv), 
= diy{VD + Z£ - Dv), 
Vt + V ■ \/v = — Vtt + iiAv — D£, 
divv = 

with initial data Z{x,0) = no{x) + po{x) and v{x,0) = vo{x). 

To obtain the desired energy estimates, we introduce new error variable = 
+p^, by the Poisson equation (2.3), we have 

= ^ , / = ^ • (2-6) 

Thus the error system can be reduced to the following equivalent system 

z^ = div{Vz^ + DE^) - A^div(£div^^ + E'^divf ) - A2div(£div£) 

- d\v{z^v^ + z\) - div(Zi;^) - A2div(^^div£'^), (2.7) 
A^[atdivi?^ - div(Vdiv^^)] + div(Z£;^) 

= -A2(9tdiv£: - Adivf) - div(5^£) - div(5^i?^) + diy{Dv^) 

- A2div(w^divf + t;div5) - A2div(u^div^^ + udivE^^), (2.8) 
+ ■ Vv^ + V ■ Vv^ ■Vv + Vtt^ - ^Aw^ 

= A^E^^divf + X^SdiYE^ + A^fdivf - DE^ + X^E^divE^, (2.9) 
divu^ = 0. (2.10) 

For the sake of notional simplicity, we set = {z^,E^,v^) and define the 
following A-weighted Sobolev's norm 

lllw^llp ^ \\{z\XE\v)\\l. + \\{z^,XE^,v^)\\lr + \\EYm- (2.11) 

The following basic inequality can be derived from Sobolev's embedding theorem 
and will be used frequently in this paper. 

Lemma 2.1. For f,g G H^{T^), we have 

\\f9\\L^<\\f\\L^-\\9\\L^<K\\f\\H^-\\9\\H^- (2.12) 
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2.1. Low order estimates. In this subsection, we derive the low order energy 
estimates from the error system (2.7)-(2.10). The first estimate is the L^{L^) 
norm of {z^,v^,E^). 

Lemma 2.2. Under the assumptions of Theorem 1.1, we have 

\\A? + \\A? + >?\\E^\? 

Jo 

< K{\\S^\\^ + + X''\\E^\\^){t = 0) 

+ K j {\\z^\\^ + + \\\-w^\\\'^){s)ds + K\'^. (2.13) 

Jo 

Proof. Multiplying (2.7) by and integrating the resulting equation over with 
respect to x, we get 

= - j DE^Vz^dx + >^^ j fdivfV^^rfa; + j vz^Vz^dx 
+ j {E^db/S + £dWE^)S7z^dx + j Zd^S7z^dx 

+ j zH^Vz^dx + \^ j diyE^E^Vz^dx. (2.14) 

We estimate the terms on the right-hand side of (2.14). For the first five terms, by 
Cauchy-Schwartz's inequality and using the regularity of D, £, v and Z, which can 
be bounded by 

e||Vz^||2 + K,\\iE^,v^, z^W + K,X^\\{E^,divE^)\\^ + K,X^. (2.15) 

For the sixth nonlinear term, by Cauchy-Schwartz's inequality and Sobolev's em- 
bedding H'^{T^) ^ L°°{T^), we get 

j z^v^Vz^dx < e||Vz^|p -hX.llw^z^lp 

<e||Vz^||2 + X,||^^|||^||5^||2 
<e\\Vz^\\^ + K,\\i^\\j,,\\zY 
<e\\VS^\\^ + K^\\\w^\\\^. (2.16) 
Similarly, for the last nonlinear term, we have 

j divE^E^Vz^dx < e||V5-^|p + i^.A^'H^-^div^-^H^ 

< e\\Vz^f + K,X'^\\E^\\l^\\diYE^\\'' 
<e||V5^||2+i^,A4||^^||2^,||div^^||2 
<e||V5^||2+ir.|||w^|||^ (2.17) 
Thus, putting (2.14)-(2.17) together and taking e small enough, we obtain 

^\\z^\f + ci\\Vz^f <K\\{z^,E^,v^)\f + KX^\\{E^,diYE^)\f 
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+ -ft:|||W-^|||'' + i^A''. (2.18) 

Multiplying (2.8) by — and integrating the resulting equation over T"^ with 
respect to x, we get 

^-^||^^||2+A2||div^^||2+ / Z\E^\^dx 

= -A^ j {dt£ - A£)E^dx - j Ez^E^dx - \^ j 6xv£v^E^dx 
-X^ j divSvE^dx + J Dv^E^dx - \^ j vE^dWE^dx 

-X^ j v^E^divE^dx - j z^E^E^dx. (2.19) 

For the first six terms on the right-hand side of (2.19), by Cauchy-Schwartz's in- 
equality and using the regularity of £, v and D, which can be bounded by 

e\\E^\\^ + KMS\i^W + K,X''\\{i\divE^)\\^ + K,X\ (2.20) 

For the seventh nonlinear term, by Cauchy-Schwartz's inequality and Sobolev's 
embedding H^{T^) ^ L°°{T^), we get 



-■^E^divE^dx < e\\E^\\^ + K.X'^Wi^diYE^W'^ 

< e\\E^\\^ + KA^\\i^\\l^\\diYE^\\^ 

pA||2 I 7^ \4||--,A||2 IIJ- ,t'A||2 



<e||£;^||2 + i^,A4||z5^||2,,||div^^ 
< fl|^^|l2 + i^,A2|||w^|||4. (2.21) 
Similarly, for the last nonlinear term, we have 

z^E^E^dx < e||E^|P + KA\z^E^\f 



<e||^^||^ + ife||5'||ioo||^^||^ 
r'A||2 , IIUT,A|||4 



< e\\E^\\^ + KS-^^IW ■ (2.22) 

Putting (2.19)-(2.22) together, choosing e small enough, and restricting A small 
enough, we get, by the positivity of Z, that 

A2^||E^||2 + 2A2||div£^||2+c2||£^||' 

< K\\{z^,v^)\f + K\\\w^\\f + KX'^. (2.23) 

Multiplying (2.9) by v'^ and integrating the resulting equation over with 
respect to x, by (2.10) and integration by parts, we obtain 

= - j DE^v^dx + X^ j v^E^dWSdx + X^ j Sv^dWE^dx 
+ X^ j v^SdhfSdx - ji^^- '^v)v^dx + j v^E^dWE^dx, (2.24) 
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where we have used the identities 

{v^ ■ S/i^)tj^dx = 0, / (w • \Iv^)i^dx = 0. 



/' 



We estimate the terms on the right liand side of (2.24). For the first four terms, 
by Cauchy-Schwartz's inequahty and using the regularity of D and £, which can 
be bounded by 

K\\{v^,E^)\\^ +K\^\\{v^,E^, div^^ )\\^+K\'^. (2.25) 
The fifth nonlinear term can be treated as follows 



[v^ ■ Vv)v^dx < K\\Vv\\loo\\v^\Y < KWv'^W'. (2.26) 

For the last nonlinear term, by Cauchy-Schwartz's inequality and Sobolev's embed- 
ding H'^{T^) ^ L°°(T3), we get 

v^E^AivE^dx < + iA*||E^divE^||2 

<^ll^5'lP + ^A4||^'lli~l|ciiv^^ll^ 

<^ll^^'lP + ^ll|w'lir- (2.27) 



Thus, putting (2.24)-(2.27) together, we get 
''-\\v'^f + 2^i\\Wv>^f <K\\{i\E^)\\^ 



dr 

+ K\^\\{v^,E^,^ivE^)f + lllw^lll^ + KX'^. (2.28) 
Combining (2.18) and (2.23) with (2.28), and restricting A small enough, we get 

+ 2;u52||Vi;^||' + (2A2 - K{\^52 + A^(5i)) ||div^^||2 
+ (c2 - K{5^ + ^2) - ^(A'<52 + A^Ji)) ll^^ll^ 
< K^\\{z^,v^)\\'' + jri|||w^|||4 + K^X' (2.29) 
for some and 82 sufficient small, which gives the inequality (2.13). □ 

Next, we estimate the L^{L'^) norm of {zt^, Vt^, Et^) by using the system (2.7)- 
(2.10). 

Lemma 2.3. Under the assumptions of Theorem 1.1, we have 

\\zi\?+\\vi\?+x'\\m' 

+ r (llVz.^lP + ||V«,^|p + \\E^f + \^\\diYElf){s)ds 
Jo 

<K{\\i^\\^ + \\i^\\^ + X^\\E^\\^){t = 0) 

+ K [ {\\{z\z^,v\v^,Vv\E\diyE^)\\'){s)ds 
Jo 
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+ K r{|||w^|||* + |||w^||n|£j^||2}(s)ds + i^A^ (2.30) 
Jo 

Proof. Differentiating (2.7) with respect to t, multiplying the resulting equation by 
z^ and integrating it over with respect to x, we get 



= - / DE^Wz^dx + y / dt{£AiY£)V~ztdx 



j dt{Zv 

+ \^ j dt{E^d\YE^)Vz^dx + j dt{z^v^)Vz^dx. 



i 



dtiz^v)Vztdx 
^'^)Vz^dx 



(2.31) 



We estimate the terms on the right-hand side of (2.31). For the first five terms, by 
Cauchy-Schwartz's inequality and using the regularity of D, £, v and Z, which can 
be bounded by 



K4\E^\\-^ + KMv^vi. 



+ K,\\{z\z^)f + KX\\{E\ Et , dWE^ , dWE^)\ \' + K,X 



A £,A 



pAm|2 



(2.32) 



For the last two nonlinear terms, by Cauchy-Schwartz's inequality, Sobolev's em- 
bedding if2(T3) ^ L°°(T3), and the inequality (2.12), we get 



/ dt{E^divE^)V5^dx+ I dt{z^v^)Vz^dx 



^A\V72A 



'sA~,A\v7;A, 



:A| 



<e||V 
<e||V5,^| 

<e\\VS^\ 



2 + K,X^\\dt{E^divE^)\\^ + if,Pt(5V)||2 
' + KX^iWE^dwEY + ll^'llioclldivE.^in 



rA||2 



+ i^,A^(||i?,-^||^.||divi^'^||^^. + ||£;-||^.||div£;i^in 



5A||2 ||~A||2 , ||;A||2 il--.A||2 
\ \h^\\V llffi 



<e\\\7~z^ 

Putting (2.31)-(2.33) together and taking e small enough, we get 
1 d 



^.lllw'lll'- 



^^\\z^\\^+cs\\Vz^\f 



(2.33) 



<K\\E^\f + K\\{v\v^)\f + K\\{z\z^)\f 



(2.34) 



+ KX'^\\{E^,E^,divE^,divE^)\f + K\\\w^\\f + KX^. 

Differentiating (2.8) with respect to t, multiplying the resulting equation by — 
and integrating it over with respect to x, we get 

yd_ 
2 dt' 



:\\E^\f + X'^WdivE^f + J Z\E^\Hx 



-I 



ZtE^EUx - X^ 



dt{dlYEv^)E^dx 



dt{dt£ - A£)Etdx - I dt{£2^)Etdx 



dt{vdiv£)E^dx + 



j dt{Dv' 



)E^dx 



dx. 
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j dt{vdb/E^)E^dx-\^ j dt{v^diyE^)diyE^dx - J dt{i^E^)E^' 

(2.35) 

For the first seven terms on the right-hand side of (2.35), by Cauchy-Schwartz's 
inequality and using the regularity of £, v, D and Z, which can be bounded by 

e||^,^||2+i^,||^^||2 + i^,||(5\sA)||2 + ;^^||(^A^^A)||2 

+ KA''\\{v\v^)\f + KA%diyE\diyE^)\f + K^- (2.36) 

For the last two nonlinear terms on the right-hand side of (2.35) by Cauchy- 
Schwartz's inequality, Sobolev's embedding iJ^(T'^) > L°°(T'^), and the inequality 
(2.12), they can be estimated as follows 

j dtii^divE^)E^dx^ J dt{~z^E^)E^dx 

<e||^,^||2+X,A4||ai(«^div^^)||2+if,||at(z^^^)||2 
<e||^,^|P+i^,A4(||5,^||?,.||div^^||^. + ||z)^||i^||div^,^in 

+ KMih\U\E^\\m+\\A\l^\\Eh?) 
<e||^,^lP+i^,A^||C,^||l,.||div^^||2,.+||^^|||.||div£:,^in 
+ if,(||z,^||^.||^^||l,.+||5^||^.||S,^in 

< ewm'' + ^.(iiiw'iir + iiiw^iini^*'in. (2.3?) 

Putting (2.35)-(2.37) together, using the positivity of Z, and taking e small enough, 
we get 

\'^^\\E^f + 2X-\\diyElf + c,\\E^f 

< i^||(f^5^^;^^;^^^)|P + ifA4i|(div^\div4^)||2 + i^A4||(S\t;,^)||2 

+ K{\\\w^\\\'^ + \\\-w^\\\^\\E^f) + K\^. (2.38) 

Differentiating (2.9) with respect to t, multiplying the resulting equation by Vf, 
integrating it over with respect to x and using divD^ = 0, we get 



■■ - j dt{DE^)v^dx + \^ j dt{£dw£)v^dx + \^ j dt{£dwE^)v^dx 
+ A^ y dt{Edbf£)v^dx - j dt{v ■ Vv^)v^dx - j dt{v^ ■ ^v)v^dx 

dt{v^ ■ Vv^)i^dx + \^ [ dtiE^divE^)ij^dx. (2.39) 



We estimate the terms on the right-hand side of (2.39). By Cauchy-Schwartz's 
inequality and using the regularity of D and £, we get 

- j dt{DE^)ildx + \^ j dti£div£)v^dx 

+ \^ j dt{£dwE^)v^dx + \^ j dt{Edw£)v^dx 

< K{\\v^\f + \\E^\\' + 11^,^11^) + i^A^||(div^\div^,^)||2 
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+ K\^\\{E^,E^)f + K\'^. (2.40) 
Now we deal with the trihnear terms involving v^jv}', v, and Vt- Using the identities 

j {v^ ■ Vv^)v^dx = 0, j{v- Vv^)v^dx = 0, 

we have 

- J dt{v^ ■ Vv^)v^dx - j dt{v ■ \/v^)d^dx - J dt{v^ ■ Vv)v^dx 

= - J ■ Vv^)v^dx - ji'^f Vv^)vtdx - j {v^ ■ Vv)v^dx - j i'"^ ' '^Vt)vtdx. 

(2.41) 

By Cauchy-Schwartz's inequality, using the regularity oft; and the inequality (2.12), 

we get 



/ 



l,,..,,o K. 



I 



<l\\vX + m^'\\\\ (2.42) 



{v,.Vv^)v^dx<^\\v^\f + ^\\vfVv^\\' 



- / {v^ . Vv)v^dx < K\\Vv\\l^\\v^\\' < Kllv^W", (2.44) 



/ 
/ 



<l\\v^\\'+K\\yv'\\\ (2.43) 



1 



2 

The last nonlinear term can be treated similarly as (2.33) 



<^\\d^\\^+K\\v^\\''. (2.45) 



J dt{E^diyE^)v^dx <^\\v^\f + ^X'^WdtiE^divE 



Am|2 



<l\\i^\\' + ^A^dl^'div^^ir + ||^^||ioo||div^,^in 
<^ll^*'ll' + lxM\E^\\U\divE^\\%. + \\E^\\%.\\divE^f) 

<^\\v^\\^ + K\\\w^\\\\ (2.46) 



Putting (2.39)-(2.46) together, we have 
':\\v^\f+f,\\Vv^f 



dr 

<K\\{v^,yv^,vi,E^,E^)\\^ + K\%{E^,E^, dvfE^ , div^i^ ) 1 1 ^ 
+ i<:|||w^|||^ + /<rA^ (2.47) 
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Combining (2.34), (2.38) and (2.47), and restricting A small enough, we get 
d_ 



-A5^\\zh?+\^\\E^\? + 5M\?) 



+ 2<S3 C3 1 1 Vz,^ 1 1 2 + /i54 1 1 Vf;,^ 1 1 2 + C5 1 1 1 1 2 + C6 1 1 divE,^ 1 1 2 
< i^2||(*\ Vv^.vti^. E\divE^)\\^ 
+ i^2(|||w^|||^ + |||w^|ini-Ei'in+if2A^ (2.48) 
for some ^3 and 84 sufficient small, which gives the inequality (2.30). □ 

Using Lemma 2.3, we can obtain the L^{Ll) norm of (V5^, Vv^, E^, Adiv.E^). 
Lemma 2.4. Under the assumptions of Theorem 1.1, we have 
||(V5^,Vu^,^^)|P + A2||div^^||2 

< K\\{z^,z^,v^,v^)\\^ +K\^\\E^\\'' + i^lllw^lll^ +ii:A^ (2.49) 
Proof. It follows form (2.29) and Cauchy-Schwartz's inequality that 

ci5^\\yA? + M^2||Vi;^||2 ^ ^2A2 - K{\^52 + A^^i)) ||div^^||2 
+ (c2 - K{5i + 82) - K{\^52 + A^Ji)) 

<-|(5i||^^ll' + '^2||i;^||^ + A^||^^||^) 
+ i^i||(5\5^)||2 + i^i|||w^|||* + /^iA4 
<K\\{z^,z^,v^,v^)\\^ +K\^\\{E^,E^)\\'' + K\\\v^^^^ 
which gives (2.49) by using Lemma 2.3. □ 

2.2. High order estimates. In this subsection we will establish the I/^(L^) of 
the higher order spatial derivatives (A^'^, Av"^, divE''^, AVdivE^). 

Lemma 2.5. Under the assumptions of Theorem 1.1, we have 

\\Az^\f + W^v^lf + ||divE^||2 + A^llVdivE^II^ 

<K\\{z^, Vz^, Vz^, v^, Vv^, Vv^, E^)\f 

+ KX^ I \divE^ \f + K\\\w^\\f + KX^. (2.50) 

Proof. Multiplying (2.7) by — A^"^, integrating the resulting equation over with 
respect to x, we get 

1|||V5^||2 + ||A5A||2 

= - j div{DE^)Az^dx + \^ j div(i?^div^ + £dwE^)/^z^dx 

+ A^y div(f divf )A5^da; + j div{Zv^)Az^dx + j div{vz^)Az^dx 

+ j dvf{z^v^)Az^dx + y div(£;^divE^)A5^da;. (2.51) 

We estimate the terms on the right-hand side of (2.51). By Cauchy-Schwartz's 
inequality and using the regularity of D,£,Z and v, the first two terms can be 
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bounded by 

e||Az^||2 +X,||(iJ^,divi?^)||2 + X,A4||(£;\divi?^,Vdiv£;^)||2 (2.52) 
and the third, fourth and fifth terms can be bounded by 

e\\AzY + K,\\i\7S\d^)\f+KA\ (2.53) 

where we use the facts that divv'^ = and divf = 0. For the last two nonhnear 
terms, using the facts that 

div(z^'D^) = Vi^f/, div(i?div£;^) = ^^Vdivi?^ + (divi?^)^, 

Cauchy-Schwartz's inequahty, Sobolev's embedding iJ^(T^) ^ L°°[T^) and the 
inequality (2.12), we have 

< e\\A2^\f + K,\\Vz^v^f + K,X'^\\diY{E^diYE^)\f 

+ if,A4(||^^||2^||VdivS^||2 + ||div^^||^i||div^^||^0 

< e\\Az^\\^ + K^WVS^W'WvYm 

+ K,X\\\E^ 11^. 1 1 Vdiv£;^ 1 12 + I |div^^ 1 11,1 1 |div^^ 1 11,1 ) 

< eWAS^W^ + K,\\\w^\\\\ (2.54) 
Putting (2.51)-(2.54) together and choosing e small enough, we have 

l\\Vz'\\^ + cr\\Az'\f 

<K\\{E^, diYE^ )\f +KX'^\\{E^, diYE^ , Vdiv^^ ) 1 1 ^ 
+ K\\{Vz^,v^)\\^ + K\\\w^\\\^ + K\^. (2.55) 

Multiplying (2.8) by divS"^ and integrating the resulting equation over with 
respect to a;, we get 

^l||div^^||2 + A2||Vdiv^^||2+ [ ZldivE^pdx 

= - j VZE^diYE^dx j div(5tf - A£:)div^^dx - j dxY{£z^)diYE^dx 

- £i^Vdiv£div^^rfa; - j vVd\Y SdxYE^dx + j v^VDdxYE^dx 

j vydbfE^dbfE^dx - \^ j v^VdbfE^dbfE^dx - j diY{z^E^)diYE^dx, 

(2.56) 

where we have used divz;'^ = and divu = 0. By Cauchy-Schwartz's inequality and 
using the regularity of £, v, D and Z, the first seven terms on the right hand side 
of (2.56) can be bounded by 

e||div^^||2 + i<:,||(^^,5^,V5^,i;^)||^ + /<:eA*||({;^,Vdiv^^)||^ + /<:eA^ (2.57) 
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By Cauchy-Schwartz's inequality, Sobolev's embedding H^{T^) ^ L°°{T^), and 
using the inequality (2.12), the last two nonlinear terms on the right hand side of 
(2.56) can be treated as follows 

- f;^Vdiv^^div.B^rfa; - j diy{z^E^)divE^dx 

<e\\divE^\f + K,X'^\\v^VdivE^\\^ + K,\\diY{z^E^)\f 
<e\ \divE^ f + K.X^Wv^Wlool \VdiYE^ f 

+ if,(||5^||ioc||div^^||2 + ||V5^||l,.||^^|i|0 

<e| |div^^ f + K,X^\\v^\\jj2 1 1 Vdiv^^ I f 

+ /^e(||5'|||,.||div^"||' + ||V5"||l,i||^'||l,0 
<e||div^^||2 + K,{1 + a2)|||w^|||*. (2.58) 

Putting (2.56)-(2.58) together and using the positivity of Z, we get 

A^^lldivS^II^ + 2A2||Vdiv^^||2 + C8||div^^||2 

< K\\{S\Vz\E\i^)\\^ + KX^\\{v^,VdivE^)\\^ 
+ K{l + X'')\\\^^\\f + KX\ (2.59) 

Multiplying (2.9) by — Aw^ and integrating the resulting equation over with 
respect to x, by (2.10) and integrating it by parts, we have 



dx 



= Vv^)Av^dx + j{v^ ■ Vv)Av^dx - X^ j E^dxwSAv^ 
-X^ j SdiyE^Av^dx - X^ j EdivEAi^dx + j DE^Ai^dx 
-X^ j E^dxwE^Av^dx + j{v^ ■ Vv^)Av^dx (2.60) 

We estimate the terms on the right-hand side of (2.60). By Cauchy-Schwartz's 
inequality and using the regularity of u, £^ and D, the first six terms can be bounded 
by 

e||Af)^||2 + K,\\{v^, Vv^, E^)]^ + K,X''\\{E^, divE^)||2 + i^,A^ (2.61) 

By Cauchy-Schwartz's inequality, Sobolev's embedding Jf^(T^) ^ L°°{T^), and 
using the inequality (2.12), the last two nonlinear terms can be treated as follows 

- A^ y E^diyE^Av^dx + j {v^ ■ Vv^)Av^dx 

< e\\Av^\f + K,\\v^ -Vv^lf + K.X'^WE^divE^W'^ 
<e||Ai;^|p + i^e||«^||ioo||Vi;^||2 + i^eA^||^^||ioo||div^^||2 

<e||Ai;^|p + i^e||«'||l,.||Vt;^||' + /^.A*||^'||l,.||div^"||' 
<e\\Av^f + K,\\\w^\\f. (2.62) 
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Putting (2.60)-(2.62) together and choosing e smah enough, we have 

+ KX'^\\{E^,divE^)\\'' + i^lllw^lll^ + KX'^. (2.63) 
Combining (2.55), (2.59) and (2.63), and restricting A is small, we get 

|((55||V5^||2 + A^lldiv^^ll^ + SeWVv^f) 

+ 55Cs\\Az^\f + (2A2 - KX% - KX'')\\ydiyE^\f + S6t^cs\\Av^\f 

+ (C8 - KS^iX" + 1) - K6eX^)\\diyE^W' 

<ii:3||(^^,^^,Vi^,t;\Vt;^)||2 + i<:3|||w^|||4 + i^3A^ (2.64) 

for some 65 and Sq sufficient small, which gives the inequality (2.50). □ 

In order to close the estimates on the right-hand side of (2.50), we need to obtain 
the uniform bounds of the time derivatives (Vz^, Vw^, AdivE'^), which is given by 
the next lemma. 

Lemma 2.6. Under the assumptions of Theorem 1.1, we have 
\\Vzi\? + X^myEh\^ + \\Vvi\? 



+ f {\\A~z^\\^ + WAv^W^ + WdiyE^W" + X^\\Vd:,yE^\\^){s)dx 
Jo 

< K{\\Vz^\f + A^lldiv^^^ll^ + \\Vv^W')it = 0) 

+ K f {\\{z\ z^ V~z\ V~z^W + \\{v\vt, Vv\Vd^W){s)ds 
Jo 

+ K [ {\\{E^,E{',diyE^,VdiyE^)\\'^){s)ds + K [ |||w^|||*(s)ds 
Jo Jo 

+ K + ||V5,^||1,. + WE^Wl. + \\diyE^\\^)}{s)ds 

+ A^ £ { 1 1 1 1 p (1 1^,^ I II,. + 1 |div^,^| II,. + 1 1 Vdiv^,^| } (s)ds + ^^A^ 

(2.65) 

Proof. Differentiating (2.7) with respect to t, multiplying the resulting equation by 
— Ai^ and integrating it over with respect to x, we get 

l|||V.?||^ + ||i.?|P 

= j |-div(D4^)+A2at[div(£:div£;-^ + £;-^div£:)]+A2at[div(£:div£:)] 
+ dtdiy{z\) + dt[diy{Zv^)]^Az^dx 

+ j {at[div(^^f;^)] + A29t[div(£;^div^^)]}A5t^rfa;. (2.66) 

We estimate the terms on the right-hand side of (2.66). By Cauchy-Schwartz's 
inequality and using the regularity of D, £, v and Z, the first integral can be bounded 
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by 

+ K,X%E^,E^,divE^,diyE^,VdivE^,VdivE^)\\'^ + K,\'^, (2.67) 

where we have used the facts divv^ = and divv = 0. For the second integral, 
by Cauchy-Schwartz's inequality, Sobolcv's embedding H^{T^) ^ L°°(T^), the 
inequality (2.12), and divi;'^ = 0, wc have 

j {dt[diy{z^v^)] + X'^dt[diy{E^divE^)]^Az^dx 

< e| I A5,^ 1 12 + [div(0 V)] 1 12 + i^.A^ I l^t [div(^Miv^^)] 1 12 
<e||A5,^||2 + i^e(||V5^li|x||i5t'||?,i + ||V5,^in|t5"ilioc) 

+ i^.A4(||^,^||io.||Vdiv^^||2 + 2||div^^||l,,||div^,^||l,i + ||^^||ioo||Vdiv^,^||^^ 

< e\\Az^\\' + K,{\\W~z^\\U\it\\l^ + W'^z^nAln^) 

+ ir,A4(||^i^||l,.||Vdiv^^||2 + 2||divE^|||,.||div^,^||^.+||^^||l,.||VdivE^^ 

< e||A5,^||2 +i^,|||w^|||^ + ;^,A^|||w^|||2(||^A||2^^ ^ ||div^,^||2,, + HVdivE^^lP). 

(2.68) 

Thus, by putting (2.66)-(2.68) together and taking e to be small enough, we obtain 
|||Vz,^|p + cio||A5,^||2 

<K\\{E^,dWE^)\\'+KM^S\Vz^,v\v^)\f 
+ K\^\\{E^,E^, divE^, dWE^, Vdivi?^, VdivE^)\\'^ 

+ i^|||w^|||4 + i^A2|||w^|||2(||4^||^, + ||div^,^||^i + ||Vdiv4^||2) +KA4. 

(2.69) 

Differentiating (2.8) with respect to t, multiplying the resulting equation by 
div£^^ and integrating it over with respect to x, we get 

^^WdiYE^lf + X''\\WdiYE^\f + [ Z\diYE^\^dx 

= - j {dt{VZE^) + ZtdivE^)diyE{'dx - X^ j dt[dtdiy£ - Adivf Jdiv^^dx 

- ^ dt[diy{z^£)]diyE^dx - j at[div(t;^divf + udivf )]div^i^da; 
+ j dt[dw{Di^)]dwE^dx-X^ j at[div(z;div^^)]div4^da; 

- j dt[diy{z^E^)]dxvE^dx - X^ j at[div({;^div^^)]div£^t^rfa; (2.70) 

We estimate each term on the right-hand side of (2.70). Noticing divu'^ = and 

divy = 0, by Cauchy-Schwartz's inequality, using the regularity of Z,£,v and D, 
the first six terms can be bounded by 

6||div4^|p + if,||(^\4\divS^)|p + X,||(z\z^Vz\V5,^)|P 

+ K,\\{d\d^)\\' + K,X'\\{i\i^)\\' + KX\\{VdwE\ ydWE^)\\' + K,X\ 

(2.71) 



QUASINEUTRAL LIMIT OF THE ELECTRO-DIFFUSION MODEL 



17 



For the last two nonlinear terms, by Cauchy-Schwartz's inequality, Sobolcv's em- 
bedding i/2(T3) ^ L°°{T^), the inequality (2.12) and divv^ = 0, they can be 
estimated as follows 

- j dt[diy{z^E^)]diyE^dx - j dt[dhf{v^dhfE^)]diyE^dx 
< e\\d:iYE^f + K,\\dt[dw{z^E^)]\\'' + K,X'^\\dt[dw{v^dwE^)]\\^ 

<e\\d:,vE^\\^ + K,{\\~z^\\l^\\diYm^^ 

+ ||V5,^||^.||^^||?,0+i^eA4(ll«t'lli=^l|Vdiv^^||2 + ||i;^||i^||Vdiv^,^in 

+ \\y~z^\\U\E^\\m)+KA\\mU\^d\^^^^ 

<e||div^,^|p + if,|||w^|||2(||div£;,^||2 + ||z,^||2^, + ||^,^||2^, + ||^^^^^ 

+ if,A2|||w^|||2(||5,^||2^, + ||Vdiv^,^in (2.72) 

Putting (2.70)-(2.72) together, using the positivity of Z, and taking e small enough, 
we get 

A^^lldiv^i^lP + 2A2||Vdiv£;,^||2 + cii j Idiv^i^pda; 

<K\\{E\EldiYE^)\\^ + KM~^\~zlV~z\Vzi)\\^ 

+ i^||(^\^,^)|P + ifA4l|(^\z),^)|P + ifA4||(Vdiv^\Vdiv^,^)|p 

+ K\\\^^\fm^Eif + \\~z^\\l. + + W^itWm) 

+ ifA2|||w^|||2(||z,^||2^, + ||Vdiv^,^||2) (2.73) 

Differentiating (2.9) with respect to multiplying the resulting equation by 
— Ai;^, integrating it over with respect to x and using divu^ = 0, we get 

1 d 



= j dt{vVv^)Av^dx + j dt{v^ ■Vv)Av^dx + j dt{DE^)Av^dx 
-X^ j dt{£diy£)Av^dx - X^ j dt{£dWE^)Av^dx - X^ j dt{Ediy£)Av^dx 

+ j dt{v^ ■ '^v^)Av^dx -X^ j dt{E^dWE^)Ad^dx. (2.74) 

By the Cauchy-Schwartz's inequality and using the regularity of v,D and £, the 
first six terms on the right-hand side of (2.74) can be bounded by 

e\\Av^\\' + K,\\{Vv\Vv^,i\v^,E\E^)\\^ 

-FXeA*||(^\^,^)|p-KXeA*||(div^\div^,^)||2 + if,A4. (2.75) 

By Cauchy-Schwartz's inequality, Sobolev's embedding H^{T^) ^ L^{T^), and 
using the inequality (2.12), the last two nonlinear terms on the right-hand side of 
(2.74) can be treated as follows 

j dt{v^ ■ yv^)Av^dx -X^ j dt{E^dxvE^)Avidx 

< e\\Avi\\^ + K,\\dt{v^ ■ Vv^)\\^ + K,X%dt{E^diYE>^)\\'' 
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+ (2A^ - KX'^ - ^7/rA*)||Vdiv^i^||^ + (ci2 - StK - KX'^{Sr + SsMdivE^^^'^ 



<e||A5,^|p + i^,(||5,^||l,.||V5^||l,. + ||z)^||i.||VS,^in 

+ i^.A^d 1^,^ 1 1^. I IdivE^ 1 1^. + 1 1^^ 1 1^. I |div^,^ I P) 
<e||Atit^||2 + l<:,|||w^|||4. (2.76) 

Thus, by putting (2.74)-(2.76) together and taking e to be small enough, we obtain 

< X| I (Vf/, VS,\ 5\ 5^ 4^) 1 12 + if A^l I (div^\ div4^)| r 
+ KX^\\{E^,E^)\\^ + KlWw^lW' + KX\ (2.77) 

Combining (2.69), (2.73) and (2.77), and restricting A is small, we get 

d_ 
dt' 

+ Ki\\{E^,E^, d\vE^ , Vdiv£;^ ) 1 1 ^ + 1 1 1 1 1 1 * 
+ K,\\\w^\\m\\z^\\],. + WVz^Wj,, + Wmi^ + WdWE^W) 
+ K4X'\\\w^\\\m\E^\\jj2 + lldivEj^ll^i + ||VdivE,^||2) +^^4A^ (2.78) 
for some 67 and 5s sufficient small, which give the inequality (2.65). □ 

3. Proof of Theorem 1.1 

In this section, we will use the energy estimates obtained in Section 2 to estab- 
lish the entropy production integration inequality and compete the proof of our 
main result. First, under the assumption of Theorem 1.1, by the standard elliptic 
regularity estimates, we have 

\\sYh^<K{\\S^\\^ + \\AS'\\% (3.1) 

\rz^rH^<K{\\S^\\' + \\A~z^\\'), (3.2) 

\v'\\l.<K{\\v'\f + \\Av^\f), (3.3) 

\v}\\%.<K{\\v^\f + \\Av^\f), (3.4) 

ll^^lll,, <if(||^^||2 + ||div^^||^3-0,s = l,2, (3.5) 

11^,^111,, <if(||^,^||2 + ||div^,^||l,._0,s = 1,2. (3.6) 

By the definitions of r^(i) and |||w^(f)||| (see the definitions (1.16) and (2.11) 

above) and using the inequalities (3.1)-(3.6), it is easy to veriiy that there exist two 
constants Ki and K2, independent of A, such that 

K^\\\w\t)\\\^ < T\t) < K2\\\^\t)\\\^. (3.7) 

Using the inequalities (2. 13), (2. 30), (2. 49), (2. 50), and (2.65), we can obtain the 
new inequality [(2.13) + 5(2.30) -l-^^ (2.49)] + j3jj(2.50) + (2.65)]. By taking 6 small 
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enough, restricting A sufficient small, and a tedious but straightforward computa- 
tion, we obtain the following relative entropy production integration inequality 

T^{t) + K [ G^{s)ds<Kf^{t = 0) + K{T^{t)f+KX'^+Kf T^{s)G^{s)ds 
Jo Jo 

+ K f {T\s) + {T\s)f}{s)ds, (3.8) 
Jo 

where G^{t) is defined by (1.17) and 

f^{t = 0) =[i|5^||2 + Wv^lf + A^ilE^II^ + \\z^f + Wv^f + A^ll^t^ll^] {t = 0) 

+ iWVz^f + A^lldiv^^^ll^ + \\yv^f){t = 0). (3.9) 

The inequality (3.8) is a generalized Gronwall's type with an extra integration 
term, we have the following result. 

Lemma 3.1. Suppose that 

f ^(i = 0) < KX^, (3.10) 

where K is a positive constant, independent of X. Then for any T G (0,Tmax), 
Tmax < +00, there exists a positive constant Aq <C 1 such that for any A < Aq the 
inequality 

r^(t) < KX^-" (3.11) 

holds for any a £ (0, 2) and < t < T . 

Since the proof of Lemma 3.1 is similar to that of Lemma 10 in [10], we omit it 
here and continue our proof of Theorem 1.1. In order to apply Lemma 3.1, we need 
to verify (3.10). In fact, by the assumptions (1.14) on the initial data (kq.Pq, ), 
we get E^{t = 0) = since the solution involved here is smooth, in particular, 
the solution and its derivatives are continuous with respect to x and t. Then, by 
using the assumption (1.14), E^{t = 0) = 0, the continuity of the solution and its 
derivatives, and the equations (2.7)-(2.10), wc get 

[\\zY + \\vY + Ws^W + WitW + A'll^t'in it = 0) 

+ (llVz.^lP + A^lldiv^^lp + \\\7v^\m = 0) < KX^ 

which gives the inequality (3.10). Thus, by Lemma 3.1, the inequality (3.11) holds. 
We easily get the estimate (1.15) by the definition of r^(f), the inequality (3.11), 
and the transform (2.6), which complete the proof of Theorem 1.1. 
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